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ADAPTED DOWNHILL SIMPLEX METHOD FOR PRICING 

CONVERTIBLE BONDS 

KATERYNA MISHCHENKO, VOLODYMYR MISHCHENKO, AND ANATOLIY MALYARENKO 

Abstract. The paper is devoted to modeling optimal exercise strategies of 
- - - the behavior of investors and issuers working with convertible bonds. This 

Q ■ implies solution of the problems of stock price modeling, payoff computa- 

O '. tion and min-max optimization. 

Stock prices (underlying asset) were modeled under the assumption of 

■ the geometric Brownian motion of their values. The Monte Carlo method 
Q ■ was used for calculating the real payoff which is the objective function. 

The min-max optimization problem was solved using the derivative-free 
Downhill Simplex method. 
' The performed numerical experiments allowed to formulate recommen- 

■ dations for the choice of appropriate size of the initial simplex in the 
pH i Downhill Simplex Method, the number of generated trajectories of under- 

[ lying asset, the size of the problem and initial trajectories of the behavior 

qh' of investors and issuers. 



O 



1 . Introduction and Problem Formulation 



1.1. Convertible Bonds. One type of securities at modern financial mar- 
ket is a convertible bond. They belong to most popular securities of modern 
financial market. This type of securities is of interest for both small devel- 
O ■ oping companies for attracting investments and for investors, since for the 

latter such bonds are highly profitable. 
X ; Strategies of the behavior of investors and issuers for early exercise deci- 

; sion must be chosen in such a way that the issuers' payofi" will be minimal 

while the investors' profit will be maximal. 

A standard convertible bond is a bond that gives the holder (investor) the 
right to exchange (convert) it into a predetermined number of stock during 



a certain, predetermined period of time (Amman et al., (2007)) 
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Convertible bonds are characterized by the following options: 
The Issuers Options 

1. Call price Kt 

This option allows the issuer to call back the convertible bond at 
the time t with the payment Kt to the investor. 

2. Call Notice Period timenotice 

Before calling back the convertible bond the issuer announces 
his intend and can call convertible bond only after the call notice 
period timenotice- During this period the investor may convert the 
bond ("force conversion"). 

The Investors Options 

1. Number of stocks n 

The investor may convert the bond into n stocks at any time during 
the predetermined period. 

2. Put price Pt 

This option allows the investor to sell the convertible bond at the 
price Pt during the predetermined period. 

Common Options of Convertible Bonds and Stocks 

1. Maturity Time T 

The expiry time of a convertible bond. 

2. Face Value of Convertible Bond N 

The predetermined price of a convertible bond which the issuer 
will pay to the investor at maturity time. 

3. Redemption Ratio k 

This is a preset percentage of the face value of a convertible bond 
which increases the price of the face value. So, the kN instead of N 
may be paid by the issuer to the investor. Usually k is equal to 1 . 

4. Price of Coupon Bond Bt 

This is the premium paid by the issuer to the investor at some 
fixed time moments during the preset period. 

5. Value of continuation Vt 

This is the price of a convertible bond at every time moment dur- 
ing the period when this convertible bond is alive. Vf is valued by 
the amount of money which the owner may get by converting the 
bond into stocks. 

6. Current Stock price Sj 

This is the current price of stocks owned by the issuer. 

7. One Option 
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The issuer or investor can perform only one action with the con- 
vertible bond, i.e. if the investor converts the bond the latter expires, 
as well as if the issuer wants to call the bond this cannot be stopped. 

Finally, the initial conditions at the time r = are: 

1. Kt>N 

The call price of a convertible bond is greater than its face value. 

2. < 

The put price of a convertible is less than its face value. 

3. n = N/{Sq • rj), where rj > I. 

n is the number of stocks obtained by conversion of a bond. This 
value is calculated under assumption that the stock price at the initial 
moment is greater than the real one. 

1.2. Problem Formulation. Objective function is as a payoff gained by 
the investor. This means that the payoff is a nonnegative number. It is 
clear that the investor wants to maximize the payoff, while the the issuer 
wants to minimize it. This payoff is based on the behaviors of the investor 
and the issuer and on current stock price at the time t. 
The problem under consideration is to find such strategies for the investor 
Convt and the issuer CalU which maximize the investor's payoff and min- 
imize the payoff payed by the issuer to investor, simultaneously. In other 
words, we have a min-max optimization problem for computing the is- 
suer's and investor's strategies. 

( 1 ) max min payofF(Co«Vf , CalU) 

Conv, Call, 

The natural choice of the boundary conditions for the investor's and is- 
suer's strategies is the following: 

(2) Convt > N 

(3) Kt > Callt > N 

We consider ([3]) only under the condition T - t > timenotice because 
otherwise the investor will not have enough time to realize his right to call 
back a convertible bond. 

Additionally, we introduce some initial settings, where we assume that we 
have a zero coupon convertible bond (no coupon payment during the preset 
period is done: Bt = 0). 
Also, we use the following preset parameters: 
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N = constant 



K 



(Kt = constant) and (Kt > N) 



So = constant. (Sq ■ n < N) 



T] > I, say 7] = 1.1 



n^N/(So-r]) 



Initial Guess: Convt - Kt + s; Calif = N + s 



time notice - const 



The maximization problem ([T]) - (j4]) will be solved by the global Down- 
hill Simplex method, see dNelder and Mead, (1964)D and stock prices are 
modeled by Monte Carlo simulations presented in ( [Garcia, (2003 )D . 

The fact that such a problem can be solved numerically using this ap- 
proach was shown in the Master Thesis ( jlsaksson, (2006)| ). 

Section 2 of this paper is devoted to the description of the methods for 
stock price generating, payoff computation, approximations of the strategy 
of the behaviors of the investor and the issuer, as well as the description of 
the Downhill Simplex method. 

In Section 3 we present the results of the numerical experiments with 
model under consideration. In these experiments we try to determine the 
best values for such input parameters as size of the simplex in Downhill 
Simplex method, the best choice for the initial trajectories, the optimal of 
trajectories used for the stock price generation and the appropriate number 
of the points for trajectories approximations. 

We finalize our work by making conclusions and giving some guidelines 
for further investigation in Section 4. 



2. Numerical Issues 

2.1. Stock Price Generating. The first step in solving the problem 
dl]) is to generate stock prices. This can be done by difi'erent methods, and 
we base our computation on the method producing the Brownian motion 
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type trajectories (see e.g. ( [Garcia, (2003)[ )). This method generates trajec- 
tories without jumps. 

The initial stock price o is given. The formula for generating the stock 
price without jumps at f + 1 time moment is: 

(4) St^i=Sf e('-^-^V2)-A+cr. ^^Z■Z.,, 

where S t is the stock price at the current time moment t\ r is the interest 
rate; 5 is the dividend yield of the issuer stocks (underlying asset); cr is the 
volatility; A = is the time step at 250 working days a year and {Zt+\,t > 
0} is a sequence of independent standard normal random variables. 

2.2. Computation of the Payoff. We use the Monte Carlo method for 
modeling a real payoff. For this purpose we generate M, a large number of 
trajectories, of the issuer stock prices according to dH), and then compute 
payoffi for each of them according to Algorithm [H 

The put price Pt is less than the face value N. In the case of maximiza- 
tion of the investor's payoff we will not take into account Pt and suppose 
that the investor does not use the possibility of choosing the put option. 

Finally, the real payoff is: 

1 ^ 

(5) payoff =j^ J] payoff 

In this study we simplify the problem by setting r - const; 6 - const 
and cr = const. 

2.3. Optimization procedure by Downhill Simplex Method. As we 

solve a min-max optimization problem ([I]), the optimization procedure is 
to be applied twice at each iteration: once as a maximizer of the objec- 
tive w.r.t. the investor's strategy and then as a minimizer w.r.t. the issuer's 
strategy. 

In other words, we find a pair of trajectories {Conv*, CalltJ which sat- 
isfy the conditions: 

(6) payoff (Conv*, Callt) = max pay off (Convt,C all t) 

Conv, 

(7) payoff(Cow*, Cfl//fH<) = min payoff (Cow*, Ca//f) 

Call, 

The optimization problem under consideration is very computational 
expensive due to the Monte Carlo simulations used of evaluation of the 
objective function, thus we do not apply the optimization methods based 
on derivatives computation. 
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Algorithm 1 Payoff (objective function) Computation 

payoff = 

flci§notice ~ 

timecheck = 

for t = : T (T = maturity • 250) do 

if 5' f > Convt then 

payoff = St (conversion or force conversion during the notice 

period) 

break 
end if 

if flagnotice = 1 then 

if timecheck = timenoUce then 

payoff = Kt (call) 

break 
else 

timecheck = timecheck + 1 
end if 
else 

if (St > Calif) and {{T - t) > timenoUce) then 

flagnotice = 1 (start call notice period) 
end if 
end if 

if t = T then 

payoff - N (the face value) 
end if 
end for 



To perform the optimization procedure in the most efficient way we use a 
derivative-free Downhill Simplex method ( |Nelder and Mead, (1964)D , which 
is also easy implementable. 

Below we give a description of this algorithm for minimizing some func- 
tion fix). 

Firstly we need to specify the initial m+ \ dimensional simplex by taking 
m points around the initial guess x^. The point with the highest function 
value fnax is called x^ax- 

The main idea of the Downhill Simplex method is to substitute a point 
with the coordinates Xmax by another point with better function value (e.g. 
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get lower f^ax in case of minimization problem). This is done by means of 
Reflection, Expansion, Contraction and Multiple Contraction. Let consider 
these functions shortly. 

Firstly, we introduce the point Xcemer (the center point of the simplex for 
current iteration): 

^ m+l 

(8) ^center ~ 7 / , ^ 

m+l ^ 

i-\ 

1. Reflection 

The point Xmax is reflected into x*j^_^ so that it lies on the opposite 
(to the point Xmax) side of the line containing Xmax and x center- The 
distance between Xcenter and x*^^^ depends on a positive constant a - 
reflection coefficient and is computed as 

(9) -^max ~ (1 ^) ' -^center ~ O! • XffiQx 

2. Expansion 

The expansion process is a prolongation of Reflection, and the 
new point x**^^ is found as: 

(10) ^max ~ y ' ^max + (1 ~ T) ' ^center 

where y{y > 1) is the expansion coefficient. 

3. Contraction 

The contraction process puts the next point between Xmax and 
Xcenter, and the new point x*^^^ can be found as: 

(11) -^max ~ P ' -^max "I" ( 1 ~ P) ' ^center 

where < < 1 is the contraction coefficient; 

4. Multiple Contraction 

Here all the points are shifted according to the following rule: 

(12) x' = (x' + x,„/„)/2, / = 1, ..m + 1 

where / is the number of points in the simplex and x,„/„ corresponds 
to the point with the minimal function value 

In Algorithm [2] we implement the Downhill Simplex method described 
in dPress etal., (1997)D . 

We denote as x' the whole strategy on /th iteration consisting of all j 
points and Xj is the j component of the strategy x'. 

So, for solving the maximization problem ^ we use the Algorithm [2l 
where we minimize the objective function -/. The problem (|7]) is solved 
by Algorithm [2l directly. 
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Algorithm 2 Downhill Simplex Method 

Choose the initial guess 

Choose the size of initial simplex k 

Choose the maximum number of iterations maxner 

Choose the termination criterion e 

for J = 2 : m + 1 do 

= x^ 
x\ - x|_j + k 

f(i) ^ ~m 

end for 

MAIN LOOP 

for i = 1 : maxiter do 

Filld fmaxt fmint fnearestmax 

if ((/ == maXiter) or {{fmax " fmiu) < e)) tlien 

break 
end if 

Reflection x„,ax to x*^^, and find /^^^ 

if fmax < fmin then 

Expansion x^^^ to x^*^^ and find the new function value /^* 
'^^fZx<fi^ax*t\itn 

max Jmax'> -^max -^max 

end if 
else 

if //Lx > fnearestmax ^^en 

Contraction X;„ax to x;;^^, and find f^^^ 

end if 
end if 

if fmax > fmax then 

Multiple Contraction 
else 

Xmax — ^max^ fmax — fmax 

end if 
end for 
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2.4. Dimension of the Problem and Investor(Issuer) Trajectory Ap- 
proximation. The trajectories Convt and Callt are the vectors of length 
250 • Maturity and 250 • Maturity - CallNoticePeriod, respectively. Since 
the lengthes of these vectors are the dimension of the optimization prob- 
lem, it is clear that such a problem cannot be solved efficiently by any opti- 
mization method. In order to reduce the dimension of the problem, solved 
by the Downhill Simplex Method we shall consider some approximations 
of the trajectories instead of the original extremely costly computable tra- 
jectories. 

So, instead of considering whole trajectories Convt and CalU in opti- 
mization procedure described above, we use a set of threshold points x,- 
the critical dates from the first possible exercise date till last exercise date 
at maturity T. The trajectories are approximated by piecewise linear func- 
tions with nodes being the threshold points x, . 

We will use 5-15 points approximation of the trajectories which means 
that the optimization problem will be 5 - 15 dimensional, too. 

Since we have a Brownian type modeling, the deviations of the prices 
from the initial state will increase while approaching the maturity. This is 
also natural for any market that the most interesting and important actions 
take place at the end. So, the distributions of the points should meet this 
requirement. According to ( |Isaksson, (2006)| ) we consider the following 



distribution of the threshold points: 

(13) = 2 

where xi = and Xm = T or x,n = T - timenotice - 1 for the investor's 
and issuer's trajectories accordingly. Also m is the number of points in 
approximation. 

Note that for the investor's trajectory the function value at the point x,„ = T 
is always equal to N. 

2.5. Short Description of the Main Algorithm. As a termination crite- 
rion we use the value of the gap. The gap is the difference between the 
investors and issuers payoff. For each iteration it is defined as: 

(14) gap = payoff(Cow*, Callt) - payoff(Co«v*, Callt ^) 

where payoff (Conv% Callt) and payoff{Conv^, Califs) are the optimizers 
for the problems ©)- ([7]) respectively. 

From the optimization point of view, due to the definition (fT4]) . the gap 
is always nonnegative. 
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The main algorithm for solution the problem ([I]) - (j4]) is sketched in the 
Algorithm m 



Algorithm 3 Body of the main loop 
gap = gapoid = 10 

while ((gap > e) and (gapou > e)) do 
gapoid = gap 

compute payoff(Conv% CalU) from ^ by one step of Algorithm [2] 
compute payoff(Conv*, CalU^) from (|7]) by one step of Algorithm [2] 
gap = payoff(Conv^, Callt) - payoff(Conv^, Callt^) 
end while 



where e is a small constant, say e - 0.0001. 

3. Numerical Results 

The numerical model described above has a quite complicate structure 
and is sensitive to the choice of the input parameters. In this section we in- 
vestigate the dependence of the performance of the method on some of pa- 
rameters. This analysis is used to validate the numerical model and choose 
the parameters which give a reasonable and fast solution. 

In each experiment we compute and compare the optimal strategies of 
the investor and the issuer. Moreover, we present the history of the behav- 
ior of the optimal value of the objective, i.e. the investor's payoff. 

While comparing the results of the numerical experiments with different 
initial parameters it is necessary to assume that due to the specificity of the 
given optimization problem the following conditions are to be fulfilled: 

• The strategies of the issuer and the investor must change their be- 
havior near the maturity time, see Section I2.4t 

• The objective function must produce dumping oscillations due to the 
nature of the min - max problem. 

Below we present the experiments where 4 parameters (initial condi- 
tions) of the optimization problem are varied. These parameters are: the 
number of generated trajectories, the size of the simplex, the initial guess 
and the number of points for approximation strategies of the investor and 
the issuer (the size of the problem). 

In each experiment the optimization problem is solved 3 times for 3 ex- 
perimental values for each of the four parameters, three other parameters 
being fixed (are from the basic set). 
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The basic set of the parameters is the following: 

(15) number of generated trajectories M = 525 

(16) size of simplex k = 3 

(17) size of the problem m = 10 

(18) initial guess e = 5 

The rest of the parameters of the problem are constant values, such as: 

Initial Stock Price = $98; 

Convertible Bond Price (Face value) = $100; 

Call Price = $110; 

Call Notice Period =10 days; 

Interest Rate = 0.05 (5%); 

Dividend yield = 0.1; 

Volatility = 0.2 (20%); 

Maturity = 2 (two years). 



3.1. Experiment 1: Different number of trajectories for Stock Price 
generating. Presented here are the results for 3 different sets of generated 



trajectories: M = 50,525 and 1000. The rest of the parameters are (fT6|) 

m. 
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Figure 1 . Strategies of Investor and Issuer. Different number 
of trajectories 
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Figure 2. Objective function history. Different number of trajectories 



Analyzing the above figures, one can see that the solution corresponding 
to the case with 50 trajectories cannot be considered to be proper as the 
number of trajectories is insufficient. Firstly, subplot 1 in Figure [I] shows 
that the behavior of the strategies close to the maturity does not change. 
This means that the amount of generated strategies has no real affect on 
the strategies of the investor and the issuer at the end of the bond lifetime. 
Secondly, Figure [21 shows that the method terminates rather fast, which is 
not appropriate for this min-max problem. 

The behavior of the investor's and the issuer's optimal strategies as well 
as the objective function history are similar in the cases with 525 and 1000 
trajectories (see Figure [2] and subplots 2-3, Figure [I])- Thus, these num- 
bers of generated stock price trajectories can be accepted for future ex- 
periments. In the basic set (fT5|) - (fT8|) we consider 525 trajectories since 
the time needed for function evaluation for this case is much shorter than 
the one for the case with 1000 trajectories, but the solution is acceptable 
thereat. 

3.2. Experiment 2: Different sizes of simplex. This experiments con- 
cerns the proper choice of the parameter k, which is the distance between 
the neighbor nodes in the initial simplex used in Downhill Simplex method. 
This experiment we run for 3 sizes of the initial simplex: k = 1,3 and 5. 
The rest of the parameters are (fT5|) and (fTTj) - (fT8)) . 
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Figure in shows that for the minimal size of simplex A: = 1 we have oscil- 
lation with small amplitude. So, the behavior of strategies of the investor 
and the issuer do not change essentially w.r.t the initial guess (see subplot 
1, Figure [3]). 

On the contrary, for maximal size of simplex k = 5, the first solution of 
the minimization has a dominating effect. In other words, the first step 
down (the solution of the minimization problem) has a big magnitude, 
which does not allow to produce sufficiently big second step up (the so- 
lution of the maximization problem). This effect manifests itself in sub- 
plots. Figure m where the functions corresponding to the strategies of the 
investor and the issuer are flat near the maturity time. 

So, we choose the size of simplex ^ = 3, which produces reasonable 
steps for both strategies (see subplot 2, Figure [3]). 

3.3. Experiment 3: Different initial trajectories. For any global opti- 
mization procedure the initial guess, which is the initial trajectories are of 
importance. We analyze 3 values of the initial trajectories e = 1, 5 and 9. 
The rest of the parameters are (fT5l)-(fT7b. 



120 F 



Results of Global Optimization 



» 110 

C 

"a> 100 

u 

°= 90 
80 

I 

110 

'I 100 
u 

°^ 90 
80 

I 

110 
■- 100 

03 
O 

i 90 
80 



* Trajectory of Investor 

• Trajectory of Issuer 
Call Price 

- Price of Convertible Bond 



50 



100 



150 



200 



250 



300 



350 



400 



450 



500 



1 1 1 1 1 1 1 1 1 


— . — . — — .. — ._. — . — . — — .. — ._. — . — . — — .. — ._. _. — ^_ 




* Trajectory of Investor 

• Trajectory of Issuer 
Call Price 




■ ■ r 1 ■ 1 ■ ■ ■ ■ r ■ ■ -1 ■ ■ ■ ■ r ■ 


- Price of Convertible Bond 





50 



100 



150 



200 



250 



300 



350 



400 



450 



500 



1 1 1 1 1 1 1 1 1 






* 






- *- 




* Trajectory of Investor 

• Trajectory of Issuer 
Call Price 






- Price of Convertible Bond 









50 



100 



150 



200 250 300 
Life time of bond 



350 



400 



450 



500 



Figure 5. Strategies of Investor and Issuer. Different initial trajectories. 

Figure [6] shows that all three experiments terminate with the same objec- 
tive function value, but give different points (strategies), see Figure [51 

The case with e = 1 requires the maximal number of iterations (almost 
twice as many as the case for e = 5 and three times as many as for e = 10). 
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The case with e = 9 seems to be not very informative, since almost 
nothing happens close to the maturity (see Figure [51 subplot 3). 

So, the most interesting cases are e = 1 and s = 5, but we choose e = 5 
in the basic set because in this case the number of iterations is twice lower 
in comparison with the case e = 1, and it produces an acceptable result. 

3.4. Experiment 4: Different sizes of the problem. Instead of using 
the whole trajectories for issuers and investors, we used the approximated 
trajectories, see Section 12.41 The amount of points in (fT3|) which gives a 
reasonable solution to the problem is the subject of investigation in this 
experiment. 

We solved our problem for different sizes of the problem: m = 5, 10 and 
15. The rest of the parameters are (fT5|) - (fT6)) and (fT8|) . 

Figure [8] shows that the objective function history for all the cases is 
almost similar. 

Nevertheless, 5-point approximation of the strategies is not sufficient. 
Since the most interesting part of the strategy is the second, it is not enough 
to have only 3 points for approximation of the second part of the strategy. 
As seen from subplot 1, Figure |7l the strategies are not smooth enough in 
the vicinity of the maturity time. 

1 5-point approximation gives very interesting results, but the size of the 
problem becomes too high as well. So, the best choice for the basic set 
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Figure 7. Strategies of Investor and Issuer. Different problem sizes 
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Figure 8. Objective function history. Different problem sizes 



is 10-points approximation of the strategies, which is the tradeoff between 
two other approximations. 
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4. Conclusions and Suggestions for Further Investigation 

In this study we considered a method for computing the strategies of 
the investor and the issuer deahng with convertible bonds. This method 
consists of two main stages: stock price generating and solution of the min- 
max optimization problem. For stock price generating we used a Monte- 
Carlo method based on the formula ©, and applied the Downhill Simplex 
method for the solution of the global nonlinear optimization problem. The 
results of our investigation allow to draw the following conclusions. 

1. The proposed method is sensitive to the number of generated tra- 
jectories of Stock Price. It means that for some small number of 
generated trajectories the method does not produce any reasonable 
solution. We suggest 500 trajectories as the minimal number re- 
quired for achieving a reasonable solution (see Section [3]T]); 



The Simplex Downhill method is sensitive to the size of the initial 
simplex. It is very important to choose the initial simplex of a proper 
size, otherwise there exists a risk to get non-acceptable solution, see 
Section [321 We recommend the size of simplex k = ?>. 



The Downhill Simplex is also sensitive to the choice of a good initial 
guess. The best choice in our experiments was s - 5 (see Section 
33]); 



4. The dimensions (size of the problem) is important in our experi- 
ments, too. Very large size of the problem requires too much com- 
putational time, but for a small size we get non-acceptable solution. 
We took 10 points (solved 10-dimensional problem), see Section 

m 

For further investigation the following is to be taken into account: 
All the results presented in this study were obtained for predetermined 
constant values such as initial stock price, call price, face value of convert- 
ible bond, etc. So, it is very interesting to run experiments with other values 
of the economic parameters. Also, such parameters as volatility, interest 
rate, dividend yield may vary during the lifetime of the bond. For example, 
they may be recalculated every day, which will make stock price genera- 
tion more complicated. Finally, the problem may be solved for nonzero 
coupon bond. 

We used Brownian type stock price generation without jumps (see Sec- 
tion [2]T]) which is one of the options. It is possible to consider some other 



stock price generation algorithms which have another nature (with jumps) 
and may give an interesting effect on the results. 



Adapted Downhill Simplex Method. 
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From the viewpoint of the optimization it would be extremely useful to 
consider another global solver, since Downhill Simplex is so sensitive to 
the choice of the initial point and the size of initial simplex. On the other 
hand, some local optimization methods may be useful, since the problem 
is constrained and the feasibility area is quit narrow. 

For more efficient strategies approximation it may be very helpful to con- 
sider other distribution of points (e.g. equidistant distribution) and other 
ways of approximation, e.g. cubic splines. 
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